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Abstract 

A relativistic equation is deduced for the bound state of two particles, by assuming a 
proper boundary condition for the propagation of the negative-energy states. It reduces to 
the (one-body)Dirac equation in the infinite limit of one of the constituent mass. It also 
has the symmetries to assure the existence of the anti-bound-state with the same mass. 
The interaction kernel(pseudo-potential) is systematically constructed by diagonalizing the 
Hamiltonian of the background field theory, by which the retardation effects are included 
r~| ! in the interaction. Its wave function at the origin(WFO) behaves properly in a manner 

sugested by the covariant field theory. 

1 Introduction 

One of the unsatisfactory nature of the Bethe-Salpeter equation for the bound state 
is that it does not reduce to the Dirac equation in the infinite limit of the one of 
the constituent mass, when the interaction is assumed to be instantaneous]]]]. We 
have to sum up all the crossed diagrams to recover the Dirac equation [f2|, which is 
impossible for the finite masses. 

Historically, the relativistic single-time equation for the two-body system pre- 
ceded the BS equation. Soon after the discovery of the Dirac equation, Breit pro- 
posed the equation of the form|| 

{H 1 (j> 1 ) + H 2 (p 2 ) + V}1> = E4>, (1) 
where H; is the Dirac Hamiltonian 



*The main part of the paper was presented in the workshop "Fundamental Problems in the 
Elementary Particle Theory" held at Nihon University in March 1995. 
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and V is a local potential. The Breit equation reduces to the Dirac equation in the 
limit mentioned above but does not have the ".E-parity symmetry", by which we 
mean that there is symmetric negative eigenvalue E for each positive one, which is 
interpreted as the bound state of the antiparticles. 

.E'-parity symmetry is a consequence of the PCT invariance, for which we assume 
the non-quantized transformations. The Dirac Hamiltonian is odd under the PCT 
transformation and the interaction Hamiltonian considered below is invariant under 
it. We therefore see that the .E-parity symmetry in the instantaneous BS equation 
is assured by a projection factor 

A ++ - A__ (2) 
in front of the potential V, which consists of the energy-projection operators 

A ^(Pi,P 2 ) = A e(Pi) A J(P 2 ) ; £,7? = + or-, (3) 

where 

Al( Pi ) = {Ei(jpi) + eH i (p i )}/2E i (p i ), 

E t ( Pt ) = (p* + m^. 

It is necessary to introduce a similar factor in any attempt at the construction of an 
improved two-body equation. 

The factor fl2|) comes from the Stiickelberg-Feynman boundary condition for the 
propagation of the negative-energy state ||. We will construct the equation for the 
unequal-mass constituents by imposing similar boundary condition and also inves- 
tigate the equal-mass equation. But, before presenting it we should restrict the 
framework of consideration. For defmiteness, we assume Abelian gauge fields in- 
teracting with the Dirac particles. We also work in the rest frame of the bound 
system(P = (E, 0)), since we are looking for the non-covariant approximation of the 
low-energy dynamics, p and x, in the following, are the relative momentum and 
coordinate, respectively, in this frame. 

2 Unequal-mass equation and its properties 

By assuming a boundary condition for the propagation of the negative- energy states, 
we deduce a new equation, which we call the Two-body Dirac equation. We first 
consider the unequal-mass constituents and assume that the mass mi is larger than 
m 2 . 
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2.1 Two-body Dirac equation 

We start with the pseudo-4-dimensional form of the equation in the momentum space 
^(p) = iSP(P,p) J V(p,q)^(g)A/(2vr) 4 , (4) 

/r)\ 

where S F (P,p) is the 2-body propagator in the lowest order and V(p, q) represents 
the interaction, which is assumed not to depend on the relative energies but is not 
necessarilly an instanteneous local potential. We impose the boundary condition that 
the negative-energy states propagates backward in time. If we retain individuality 
of the constituents and use the usual Feynman propagator the factor @ results. 
However, we can choose the other possibility in which we incorporate the idea that 
the bound two bodies should be treated as a quantum-mechanical unityfj Since m\ 
is larger than m 2 , the free part of the Hamiltonian has the same sign as that of 
the particle 1 in the CM system. Let us then modify the boundary condition as 
follows: A bound two-particles state propagates backward in time if their net energy 
is negative. By assuming it, we have the lowest propagator 



SP(P,P) - E _ pQ _ Hi{ _ p) + ^ + pQ _ HM + ie5 ^ n lhl (5) 

where Ai + A 2 = 1 and the limit 5 — > +0 is assumed. 

After integrating out the redundant degree of the freedom in (Bj), we get the 
required equation 

{Fi(-p) + H 2 (p) +$>A«,V}V> = Ei(}, (6) 

which reduces to the Dirac equation in the infinite limit of m\. It is easy to see the 
.E-parity symmetry of this equation. 

If we would apply our equation to the scattering state, we shall have, from the 
time-dependent equation, the conserved probability density 

p(x, t) = *f){x, ty £ ^ v ip(x, t), (7) 

er) 

which is in accord with the boundary condition that the negative-energy state prop- 
agates backward in time carrying the negative probability density But, it does 

*We can establish the concept of individuality in the quantum mechanics only through obser- 
vation, which brings about a subtle point to the bound system even in infrared-free theories: To 
detect an individual one in the bound constituents, we need to separate them by applying the 3rd 
interaction, which inevitably destroys the original state. So there is no reason why we have to apply 
the free propagator individually to each constituent. 

2 We note that we can interprete the causal propagator of a Dirac particle, applied to one-particle 
wave function with the negative energy, also as carrying the negative probability density 
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not necessarily provide the normalization condition for the bound state. For the 
scattering processes, the projected wave function A £V tp corresponds to the physical 
state of the (free)particles with the positive or negative energy E. And the above 
interpretation of the probability current actually says that the state with the negative 
E is carrying the negative probability density. However, for the bound state with the 

positive eigenvalue E, A ip or A_ + tp is merely a negative-energy component in the 

representation in which the energy of the free particle is diagonal. It is like a small 
component of the one-body Dirac equation. 

Taking the above consideration into account, we restore the probability interpre- 
tation of the wave function and normalize it by assuming the probability density 

p{x) = if>{x,t)ty(x,t). (8) 
Observables except for the Hamiltonian are self-adjoint under this metric: 

(0,6^) = (O0,VO- (9) 

The Hamiltonian is the operator ruling the time development of the system and 
is modified by the factor S £7? £A e?7 . Though it is not a self-adjoint operator, its 
eigenvalue is proved to be real if the inner product (|i~3|) below exists. 



2.2 Green's function and the vertex equation 

The Green's function G for (|6|) satisfies the operator equation 

{E — Hi — H 2 -J2eAe V V}G = 5>A £7) , (10) 

and 

G {E - H x - H 2 - V eA £V } = £ sA £V . (11) 

erj erj 

In the momentum representation, it can be written, by using the eigen-f unction Xnip) 
of (H), as 

G(p, p') = ~\F7p TT\Xn(p)Xn(py + continuum, (12) 

n ^n{E - M n ) 

where M n is an eigenvalue and N n is a normalization factor defined by 

(13) 

erj 

When one of the constituents(labeled with 2) is in the category of the antiparticle 
of the other, there can be an annihilation process for which the unamputated-decay- 
vertex $ is given by 

$ = C 7 G, (14) 
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where 7 is the lowest vertex and C is the charge-conjugation matrix of the particle 
2. 

$ satisfies the vertex equation 



$(£ -H x -H 2 - V^2eA ev ) = C 7 £eA, 



(15) 



and the amputated vertex is 



r = 7 + CW. 



(16) 



We can determine the renormalization constant Zi for the wave function at the 
origin (WFO) from fllBD and ([If]), if we need it J?]. 

2.3 Interaction Hamiltonian 

We have, so far, not specified the interaction Hamiltonian(quasipotential). In this 
section, we investigate it for the one-(Abelian)guage-boson exchange in the Coulomb 
guage as an example. For the instantaneous part of the interaction, V is obvious. 
For the remaining part, we can specify the quasipotential in a clear way from the 
background field theory. We have already shown, for the Salpeter equation, that 
the quasipotential from the one-boson exchange is given through the diagonaliza- 
tion of Fukuda, Sawada, Taketani[|J and Okubo[p|(FSTQ)p!0|. However, we first 
have to correct some error in Ref . HlOf -PI Namely, we employed the usual Fock space 
and reinterpreted the matrix elements of the interaction Hamiltonian including the 
negative-energy indices as the one in this space. The guiding principle was the hole 
theory. But it brings the procedure into confusion, since we have revived the nega- 
tive energy in Eq.(^j). The correct choice is to generalize the Fock space by ignoring 
the hole theory. The hole theory is partially recovered afterward by including the 
projection factor in Eq.@. 

We show only the lowest one-boson exchange potential in the following. We first 
introduce the generalized Fock subspace of the free constituents, the bases of which 
are denoted by |e, r], p), where e and r\ are the signs of the energies of the particles 1 
and 2 respectively. We then diagonalize the Hamiltonian in the Schrodinger picture 
by using the FSTO method. The second-order boson-exchange potential in this 
subspace is given by 



where q = p — p' '. The retardation effects are included in this equation. 

3 The error is only conceptual for the Salpeter equation. The result of Rcf.pOf is correct. 




2 L g 2 - {eE x {p) - e'E^Y q 2 - {r]E 2 (p) - r)'E 2 (p')} 



(17) 
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2.4 WFO of the 1 S state 

Let us next study some fundamental feature of the equation. When it is applied to 
the system in which the pair annihilation of the constituents can occur, an important 
physical quantity is the wave function at the origin(WFO). For example, the decay 
amplitude of the pseudo-scalar Qq meson via a weak boson is proportional to the 
average WFO Tr{7 5 7oV'(0)}, where ?/>(0) is the charge conjugated (with respect to the 
particle 2(g)) WFO. We investigate, in Appendix, the asymptotic behaviour of the 
momentum-space wave function by using the method given in Ref.p]]. We assume 
instantaneous exchange of a gauge boson[]. The average WFO thus obtained is finite. 
This result is consistent with consideration on the covariant field theory. We note 
that the average WFO becomes divergent in the limit of the one-body Dirac equation, 
for which we have the renormalization procedure [|TT[ . 



There are many "two-body Dirac equation" proposed. An interesting one from 



the point of view of the present paper is the one by Mandelzweig and Wallace [12 
Instead of redefining the two-body propagator, they intended to include the effects of 
the higher-order interaction (the crossed Feynman diagram) and obtained an equa- 
tion which has the proper one-body limit and the .E-parity symmetry. An important 
difference from our equation is in the average WFO considered above. It is finite in 
the Coulomb model but divergent if the transeverse part of the gauge-boson exchange 
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2.5 On the three-body equation 

The next comment to be made is on the three-body equation in which the interaction 
is assumed to be composed of two-body interactions. It is almost straight forward 
to deduce the three-body Hamiltonian. However, we have to examine how to define 
the projection factor in front of the quasipotential, since the sign of the energy of 
a interacting sub-(two body)system depends on the reference frame. This shows 
that the applicability of the single-time formalism is more restrictive in the three- 
body system than in the two-body system. The effective range of the velocity of 
the subsystem relative to the three-body CM system cannot be large. Under this 
restriction, we define the projection factor in the rest frame of the subsystem, which 
is, then, transformed into the rest frame of the three-body system. A part of the 
Hamiltonian which represent the interaction of the particle 1 and 2, for example, 
becomes 

H 12 = J2 £A l(Pi)Vi2(p v P 2 ), (18) 

£ 

when m 1 > m 2 . 

4 The analysis in the Appendix cannot be applied to the retarded interaction. 
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Brown and Ravenhall pointed out that the three-body equation of some category^ 



does not have proper solutions with the normalizable eigenfunction fll5[ . This phe- 
nomenon is called "continuum dissolution(CD)" . The criteria for CD are as follows: 
(a) The equation decouples into the independent ones when a part of the interaction 
is switched off. (b) The subequations have solutions of the positive- and the negative- 
energy continua. Our three-body equation does not fulfill (a) and is free from CD. 



3 Equal- mass equation 

So far we have considered the unequal-mass equation. In this section we investigate 
the equal-mass limit of it, for which the exchange symmetry is an issue. We first 
note that the projection factor in front of the interaction term of @ includes a part 
which violates this symmetry: It is shown that 

A (v) = A+_ - A_+ (19) 

and the Heisenberg's exchange operator 

p H = ~(i + <n ■ a 2 )(i + Pl ■ p 2 )P M (20) 

anticommute for ni\ = rri2, where the operator Pm exchanges the momenta (or 
coordinates). violates the symmetry since the remaining part of the Hamiltonian 
is commutable with P H . For equal mass, we have two equations. One is the equation 
(|6[) with mi = m 2 and another is obtained by assigning a minus sign in front of A< v > 
which is the equal-mass limit of the equation with m 2 > m\. It is the conjugate 
equation of ([]) in the sence that (^) is converted into it by the transformation Ph. 
It is easy to show that these equations have the common eigenvalue spectrum: If 
an eigenfunction \n of © belongs to some eigenvalue M n , PhXu is the solution of 
the conjugate equation with the same eigenvalue. However, \ n does not have the 
definite P#-parity. 

A way to recover the exchange symmetry is averaging the Hamiltonians of the 
equation (H) and the conjugate one. The resulting is the Salpeter equation. The 
eigenvalue of this equation is different from the above M n . The difference is, however, 
small since it is of the 4th order in the symmetry-breaking part of the Hamiltonian. 
The quarkonium phenomenology for the equal-mass constituents is given in Ref. (TO . 
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Appendix 

We examine the asymptotic (p — > oo) behavior of the momentum-space wave function 
and show that the average WFO Tr{7 5 7o^(0)}/v^ is finite.^ 

There are 4 partial amplitudes h £r) (p) in the 1 So state, with which the wave 
function is expanded as 

.W^E^WWWp)^)'^ (2D 

where c\/ 2 = — C-1/2 = l/v2 and the spinors it and v, for the particle 1 and 2 
respectively, are defined in 0. 

The average WFO for the annihilation decay through the axial- vector current is 
given by 



J2i\/( E i + srnx)(E2 + r]m 2 ) - er 1 J(E l - em l )(E 2 - 7]m 2 )}h £r) (p)p 2 dp. (22) 



ST) 



We first assume the Coulomb potential. The partial-wave equation for the 1 5'o 
state is given by 



E"rf 

x q - \ 1 w r 

p l E l (p)E 1 (q)E 2 (p)E 2 ( q y 

J Al A 2 1 I a1 a2 



{E - eE 1 (p) - v E 2 (p)}h £V (p) = -e^- ]T / dq 

e'-q' 

\{A\ el A 2 m , + se'Tjrj'Al^Al^Qoiz) 
+ {ee'A^A^, + vrfA^AL^yQ^heviq), (23) 
where z = (p 2 + q 2 )/2pq and Qe(z) is the Legendre's function. A l ££ , is defined by 

A\ £ , = y/(Ei(p) + e mi )(Ei(q) + e'rm). 

The asymptotic behavior of the wave function is determined from the integral 
region near the infinity. We then expand the both sides of (^3|) into the series of 
1/p and 1/q. We assume the power behavior of the wave function for large p. The 
independent amplitudes are chosen to be Ha{p) = J2 £ h ££ (p), h B {p) = J2e £ hes(p), 
hcip) = J2e^e-e(p), an d ho(p) = J2e £ ^e-e(p), which are expanded, in the high- 
momentum region, in power series of 1/p: 



hx{p)=Y.C n xP 



n 



'See Ref.|l(| and references therein, for the Salpeter equation. 
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Integrals on the right-hand side can be done if we neglect infrared- divergent terms 
which are irrelevant to the leading asymptotic behavior. Now, we can determine the 
asymptotic indices Px's from consistency]^]: We get, for Ka and h B , 

2C° A p-^ - EClv^*- 1 = -C\-^— cot(^)p-^ (24) 

7r I- pa 2 

and 

X%p-P° - EC\v~^ = -C° B ^-^\ tan(^ B )p^, (25) 

7r p B (2-p B ) 2 

where the terms of the higher power in 1/p are neglected. If we neglect the second 
term in the left-hand sides of (p4|), we find Pa in the range 1 < Pa < 2 Q and get 
Pb — Pa + 1 from (p5|) . We obtaine another series by neglecting the second term in 
(|25|). For this, P B is found to be in the range 2 < P < P B < 3, where the lower 
bound /?o corresponds to the upper bound 4 fir of a above which the index Pa from 
( |24|) becomes complex. /3a of the second series is given by Pa = Pb + 1- 

The asymptotic amplitudes /ic* and /ip are determined dependently on /i^ and 
/iB. We get, for the minimum indices 

Pc = mm(p A + 2,p B + 1) (26) 
Pd = Pa + 1. (27) 
We see that the average WFO ( ^2|) is finite, because 

P B > I and /3c > 2 

hold for the asymptotic amplitudes. This conclusion is valid even if the instantaneous 
exchange (transverse part) of the gauge boson is added. 
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